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Abstract 

We study a confined mixture of bosons and fermions in tlie quanta! degeneracy regime with 
attractive boson-fermion interaction. We discuss the effect that the presence of vortical states and 
the displacement of the trapping potentials may have on mixtures near collapse, and investigate 
the phase stability diagram of the K-Rb mixture in the mean-field approximation supposing in 
one case that the trapping potentials felt by bosons and fermions are shifted from each other, as 
it happens in the presence of a gravitational sag, and in another case, assuming that the Bose 
condensate sustains a vortex state. In both cases, we have obtained an analytical expression for 
the fermion effective potential when the Bose condensate is in the Thomas-Fermi regime, that can 
be used to determine the maxima of the fermionic density. We have numerically checked that the 
values one obtains for the location of these maxima using the analytical formulas remain valid up 
to the critical boson and fermion numbers, above which the mixture collapses. 

PACS numbers: 03.75.-b, 03.75.Mn, 03.75.Lni, 32.80.Pj 
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I. INTRODUCTION 



Recent experiments on degenerate Fermi-Bose mixtures have opened the possibihty of 
studying in a direct way the effect of quantum statistics in Bose-Einstein condensates 
(BECs). In practice, the direct evaporative coohng techniques used to obtain BEC are 
not apphcable in a Fermi gas, as the Pauh exclusion principle forbids s-wave collisions be- 
tween fermions. This difficulty has been overcome using a gas of bosons as a coolant, which 
has given further relevance to the study of these mixtures. 

One of the systems studied in recent experiments is the ^Li-^Li mixture 

HQ. 

It is 

characterized by having a positive boson-fermion scattering length an order of magnitude 
larger than the boson-boson one m. Due to this repulsive boson-fermion interaction, the 

n 

system does not exhibit a large overlap between the two species yj. Moreover, this mixture 
may undergo a two-component separation p|, Both facts conspire against having the 
fermions well inside the boson cloud, which is desirable to obtain an efficient sympathetic 
cooling. In contrast, this desirable overlap can be achieved in ^°K-^''Rb mixtures due to the 
large attractive boson-fermion interaction. This mixture has been recently obtained 
and it has been shown that if the particle numbers are above some critical values the system 
collapses 0]. 

From a theoretical point of view, a fairly amount of work has been devoted to the study 
of boson-fermion mixtures In particular, a systematic study of the structure of 

binary mixtures has been performed in Ref. jsi, where all possible sign combinations of 
scattering lengths between boson-boson and boson-fermion s-wave interactions have been 
discussed. The purpose of our work is, firstly, to present an analysis of the location of the 
minima of the effective potential felt by fermions submitted to a large number of condensate 
bosons when the boson-boson interaction is repulsive, and secondly, to extend this study to 
systems with a large number of fermions with an attractive boson-fermion interaction. 
An essential characteristic of superfluid systems is the occurrence of quantized vortices 
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1^|. Quantized vortices in BECs were first produced experiment allv by Matthews et al. 



and have been studied since in detail (see Refs. |12l . Il3l ll^ and Refs. therein). It 
may thus be interesting to see what is the physical appearance of such vortices, arising in 
the bosonic component, in the presence of a fermionic cloud which is in the normal (non 
superfiuid) but quantum degenerate phase. This is another aim of our work, which bears 
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some similarities with the description of quantized vortices in ^He-^He nanodroplets recently 
addressed jl5|. 

This work is organized as follows. In Sec. II we describe the mean-field model we have 
used. In Sec. Ill we present an analysis of the effective potential felt by the fermions 
when they are inside a large condensate that exhibits a boson-boson repulsive interaction, 
and the number of fermions (Np) is much smaller than the number of bosons (Nb). We 
shall consider the case in which the boson and fermion external confining potentials are 
displaced from each other, which may be caused by gravity for instance, and also when 
bosons are in a vortical state. In Sec. IV we present the results obtained by solving the 
mean-field coupled equations for different Nb and Np values up to the critical values where 
collapse of the mixture occurs. A summary is given in Sec. V. Finally, in the Appendix we 
derive some expressions using a scaling transformation which, on the one hand, constitute a 
generalization of the virial theorem, and on the other hand are especially useful for testing 
the accuracy of the numerical procedure. 

II. THE MEAN-FIELD MODEL 

We consider a mixture of a Bose condensate (B) and a degenerate Fermi gas (F) at zero 
temperature confined in an axially symmetric harmonic trap. Assuming that the minimum 
of the trapping potential felt by each species may be displaced in the z direction in a value 
di (i=B,F), the confining potentials in cylindrical coordinates are 



uj',j,r' + ujiB{z-dBy (1) 



u'^pr' + ui^iz-dpY , (2) 



Vb = \mb 

and 

Vp = -Mp 
2 

where UrB, ^zB and UrP, oJzp are the trapping radial (axial) angular frequencies for bosons 
and fermions, respectively. Mb and Mp are the corresponding masses. 

Since the number of fermions we consider in the numerical calculations is fairly large, the 
fermionic kinetic energy density can be written in the Thomas-Fermi- Weizsacker (TFW) 
approximation as a function of the local fermion density np and its gradients. For fully 
polarized spin 1/2 fermions, it reads: 

rp{r) = %^'?''nT + P^^^ (3) 
5 np 



and the fermionic kinetic energy is 



F 



2M, 



F 



J dfrp{f) . (4) 



The value of the (3 coefficient in the Weizsacker term is fixed to 1/18. This term contributes 
httle to the kinetic energy, and it is usually neglected However, its inclusion 16| has 
some advantages, see below. We refer the interested reader to Ref. [IT] for a discussion on 
the accuracy of the TFW approximation (see also Refs. therein). 

Neglecting all p-wave interactions, the energy density functional for the boson-fermion 
mixture at zero temperature has the form 

n'^ 1 

ZIVl Q Z ZIVI p 

where hb = I^E'P is the local boson atomic density. The boson-boson and boson-fermion 
coupling constants qeb and qbf are written in terms of the s-wave scattering lengths 
and asp as gsB = ash^ / Mb and qbf = clbfTi^ / Mbp, respectively. We have defined 
Mbp = 2MbMp/{Mb + Mp). 

When bosons sustain a quantized vortex line along the z-axis, the condensate wave func- 
tion can be written as \1/ = |\l/|e*™'^, where m = 1, 2, 3 ... is the circulation number, yielding 
for the kinetic energy density 

t2 1-2 +2 2 

|v^r = 777:r(V|^lf + TTir^-#. (6) 



2Mb ^Mb 2Mb ' 

We have considered only singly quantized vortices, that is m = 1. If a vortex is present in 
the condensate, bosons fiow around the vortex core with quantized circulation, which yields 
the centrifugal term in the kinetic energy. We assume that the Fermi component is not 
superfiuid, and consider that it is in a stationary state. This situation could be achieved 
experimentally waiting enough time after the generation of the vortex in the condensate, to 
let the drag force between bosons and fermions to dissipate. 

Variation of S with respect to \1/ and np keeping Nb and Np fixed yields the following 
coupled Euler-Lagrange (EL) equations 

/ h'^V^ h^m^ 1 \ 

+ ^B + 7777 — o + 9bb nB + gsF np] ^ = ^ib^ (7) 



2Mf 



2Mb 2Mb 



(6vr^)2/^nf + - 2^^np 

np 



+ Vpnp + Qbf nB np = fipnp , (8) 



where hb and j2p are the boson and fermion chemical potentials, respectively. Then, the 
ground state (m = 0) or a vortical state (m = 1) are obtained by solving the Gross-Pitaevskii 
(GP) equation for bosons [Eq. ((Tj)] coupled to the Thomas-Fermi- Weizsacker equation for 
fermions [Eq. (jHj)]. 

The inclusion of the Weizsacker term in the energy density has the major advantage that 
it yields a EL equation for np that is well behaved everywhere, avoiding the classical turn- 
ing point problem when this term is neglected (Thomas- Fermi approximation). Moreover, 
solving Eq. (jH)) is not more complicated than solving the GP equation. This can be readily 
seen writing the later in terms of n^: 



n {VuB^ 1 



2M, 



B 



ub 2 



7;^nB 



+ Vb + 7777 5" + 9bb tib + Qbf Ub rip = I^bUb , (9) 

zMb r-^ 



which is formally equivalent to Eq. (jHl). We have discretized these equations using 9- 
point formulas and have solved them on a 2D (r, z) mesh using a sufficiently large box. 
The results have been tested for different sizes of the spatial steps (we have mostly used 
Ar = Az ~ O.lyum). We have employed the imaginary time method to find the solution 
of these coupled equations written as imaginary time diffusion equations 3|- After every 
imaginary time step, the densities are normalized to the corresponding particle numbers. To 
start the iteration procedure we have used positive random numbers to build both normalized 
densities. This avoids to introduce any bias in the final results. We have also checked that 
the solutions fulfill the generalized virial theorem deduced in the Appendix. 



III. EFFECTIVE FERMION POTENTIAL 

Before we present the numerical results obtained by solving Eqs. ((HI) and 0, it is useful 
to analyze the features of the effective potential felt by a small number of fermions in the 
presence of a boson condensate, paying special attention to the location of its minima. We 
are interested in studying the cases in which either the condensate hosts a vortex line or there 
exists a displacement between the minima of the boson and fermion external potentials. The 
later situation is routinelly met in the experiments -a gravitational sag in the z direction-, 
in which case the displacement is giXjio^p — l/cj^^), being g the acceleration of gravity. 

An interesting issue is that the analytic expressions we obtain in this Sec. for the location 
of the minimum of the effective potentials remain valid, as we will show in the next Sec, 
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when the number of fermions and bosons are similar. These locations coincide with the 
positions of the maxima of the fermion density and are relevant because the collapse of the 
mixture originates precisely around them. This will be discussed in the next Sec. 

The effective fermion potential has two main terms, the external potential and the mean- 
field term arising from the interaction with the boson condensate, which is proportional 
to ub- We consider a large condensate in the Thomas- Fermi (TF) regime with positive 
scattering length and we assume that its density profile is not affected by the fermion 
presence. To obtain we may thus use Eq. ^ with np = and neglect the first kinetic 
energy term. 



A. Effect of a shift in the minimum of the external potentials 



For simplicity, we restrict ourselves to the vortex- free case m = 0. Without loss of 
generality, the origin of coordinates can be fixed at the minimum of the boson trapping 
potential, and we will assume that the displacement is only in the z direction. Thus, the 
effective trapping potential experienced by the fermions is 

1 



2 2,2/ J 

oj^p r + oj^p [z - dz 



+ Qbf riB 



(10) 



where dz = dp — ds is the z displacement between the trapping potential centers. Assuming 
that the density profile of the condensate is not affected by the interaction with fermions, the 
number of bosons is large, and that the interaction between bosons is repulsive {qeb > 0), 
one obtains from Eq. setting m = the boson density profile that corresponds to the 
uncoupled TF density 



9bb 



^ib--Mb [ujIb + ^Ib z^) 



e 



Hb-^Mb {uj^b + ^Ib z^) 



with 0(x) = 1 if X > and zero otherwise. The TF condensate boundary is given by the 
ellipsoid arg{Q) = 0. For an axially symmetric trap it yields the well-known TF radius 
of the condensate i?,- = 



'2fiB/{MB(^iB) with i = r, z in the radial and axial direction, 
respectively. Replacing hb into Eq. ()1U|) and defining the dimensionless parameter 



aBMBFMpuj'^ 



(12) 



iF 



6 



with i = r,z, it follows that the effective fermion potential inside the Bose condensate -where 
ub [Eq. (fTTj) ] is positive- is 



1 



J 2 QBE 



whereas outside the condensate 



^ - u^Fr' + oo'ipiz-d^y . (14) 



2 

This effective potential can be viewed as having a renormalized frequency inside the boson 
condensate, a feature already discussed in Ref. for a mixture with concentric external 
potentials. This model has been called the double-parabola model by Vichi et al. |s|, and it 
has also been used by Capuzzi and Hernandez However, in these works the possibility 
of a displacement between the minima of the potentials has not been considered. 
The extremum of V^-^^ inside the Bose condensate is attained at the point 

rf^ = «,<,<) = (0,0,7,-M.) . (15) 

— * 

Note that if the shift between the external potentials is ci = {d^, dy, d^)-, the effective potential 
has its extremum at d' = {■y'^d^., 'jy^dy, ■y~^dz). Depending on the signs of 7,, d' may 
correspond to a maximum, a minimum or a saddle point. 

If the boson-fermion interaction is attractive ((7^^ < 0) as for ^°K-^'^Rb mixtures, 7i > 1 
and the point d' is a minimum. Depending on the positions of d and d' with respect to 
the ellipsoid arg{Q) = 0, with the argument of taken from Eq. (jllj) . there are three 
possibilities: a) If d is inside this ellipsoid, fermions view only one minimum at d'; b) If d is 
outside and df is inside, there are two minima at d and c?; c) If both points are outside the 
ellipsoid, there is only one minimum at d. 

B. Vortices 

We consider now the effect of a quantized boson vortex line in the fermion distribution, 
without gravitational sag. In this case, the effective fermion potential is 

= Vp + qbf nl = ]^MF {uIp + ulp z^) + qbf u^b , (16) 

where is the boson density hosting a vortex line along the z axis. This density is zero at 
r = 0, and reaches its maximum value on a circle of radius Tq around that axis. In the TP 



approximation, can be derived using Eq. (171) keeping the centrifugal term proportional 
to m, and may be approximately written as [l^ 



n 



B 



9bb 



IJ'B 



h^rn? 1 
2Mb r2 



B 







h^w? 1 

2Mb 



(17) 



with ro = \Jhni/{MBiOrB)- Using we calculate the mimimum of V^^^ inside the conden- 
sate, which is located at = and 



(18) 



V MBUJrB 

Thus, the minimum is attained at a circle of radius r'. We will see that for the ^°K-^''Rb 
mixture r' is very close to the radius at which the boson density has its maximum. 



IV. NUMERICAL RESULTS 



The system under consideration is a confined ^°K-^''Rb mixture. We have assumed spher- 
ically symmetric traps for bosons and fermions, with trap frequencies ujb = 27r x 100 Hz 
and Up = \J Mb / Mp ujb- When a displacement dz is introduced between the minima of 
the external potentials, the mixture has only axial symmetry around the z axis. We have 
numerically solved Eqs. (jH)) and using the set of scattering lengths reported by Modugno 
et al. 0, namely = 98.98 Cq, and a^F = — 395ao, being the Bohr radius. The dimen- 
sionless parameters introduced in Eq. (fT^ are 7"^ = 7"^ = 0.136, and (1 — 7^"^)^^"^ = 0.964 
[Eq. mi 

We display in Fig. Q several boson and fermion density profiles as a function of z, 
considering an arbitrary displacement dz = 10 fim. They all correspond to Nb = 10^, but 
for three different fermion numbers: Np = , Np = 10^, and Np = 2.5 x 10^. The attractive 
boson-fermion interaction produces an enhancement of the density of both species in the 
overlap volume. However, this effect is reduced with respect to the concentric case due to z 
displacement 0] • 

Using the TP approximation, the radius of the Bose condensate can be estimated as 
Rb = {I'^NBCiB / ihoY^^clhO) with auo = \Jf>'/ Mbojb being the oscillator length of the 
boson trap; this yields Rb ^ 6/im. Following the analysis performed in Sec. HI. A, and 
using 7~^ = 0.136, we find that d'^ = 1.36;um < Rb; thus, within that model there are two 
minima in the effective fermion potential. One is at d = 10 /xm z, i.e., beyond the TF radius, 
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and the other is inside the boson cloud, at d! = 1.36 /im z. It is interesting to note from Fig. 
Qthat in all cases, the maxima of the fermion density are precisely at z = d'^ and z = dz, 
which we have displayed in the inset as vertical lines. 

We thus see that for Nf/Nb equal to 1% (solid line) or even 25% (dashed line), the 
maxima of the fermion density appear at the values obtained using the model of Sec. III. A. 
This can be understood by solving Eqs. ((Tj) and (jH)) with m = (i.e., no vortex), without 
neglecting the effect of fermions in the GP equation as it was done in the simplified model 
of Sec. Ill, and assuming that all the gradient terms in both Eqs. can be neglected. This 
assumption is justified for the number of bosons and fermions we have used: in the case of 
bosons, this is what the TF approximation is about, and in the case of fermions, the gradient 
term is a correction to the leading term proportional to n^J^ . Using Eq. ((7j), we obtain the 
following expression for n^: 

ub = ^— [fiB -Vb - Qbf np] . (19) 

9bb 

Replacing it into Eq. (jH)) and deriving the resulting expression with respect to z we get 

3Mf OZ QBE 9bb oz 

The extremum of the fermion density is found by setting dnp/dz = 0, with np ^ 0. After 
some algebra we find that the following equation has to be fulfilled: 

Mpulj,{-fzz-dz) = , (21) 

whose solution coincides precisely with the value z = d'^we have found in the previous Sec. 

Whereas the maximum of the fermion density is located aX z = d'^ irrespective of the 
value of the dilution (Np/Nb), it can be seen from the numerical results that the maximum 
of the boson density moves from z = towards z = d'^ as Np increases. In fact. Fig. 
^ shows that even for a rather small Np value {Np = 10^, Nb = 10^), the shape of the 
condensate differs from the parabolic-type profile yielded by the TF approximation for a 
fermion-free condensate, Eq. (fTTjl . It can be also appreciated in the inset that for a dilution 
of Np/Nb = 25%, a fairly large number of fermions remains without mixing (long tail in 
the fermion density profile outside the Bose condensate). 

We plot in Fig. |21 several density profiles for configurations hosting a singly quantized 
vortex line in the case of no displacement between the trapping potentials. We have taken 



Nb = 10^, and have considered three different fermion numbers, Np = 0, 10^, and 1.5 x 10^. 
When Np = 0, the maximum of the bosonic density is located in the z = plane, in a 
circle of Tq = 1.13 /im radius marked with a vertical line in the graph. Fermions are peaked 
around r' = 0.96 rg, as we have obtained using the simplified model of Sec. III.B with 
(1-77^)1/^ = 0.964. 

It may be seen that in the presence of a quantized vortex, the position of the maximum 
of boson and fermion densities are very close irrespective of the dilution Np/Ns- It is 
interesting to note that increasing Np at constant Nb, the fermion and boson density profiles 
become sharper due to their mutual attraction, eventually collapsing for a critical number 
of atoms. 

Following a procedure analogous to that used before, we can justify the position of the 
maximum of the fermion density in the presence of a vortex in a TF condensate. From Eq. 
((Tj) we obtain 

1 h'^m'^ 1 

ns = [f^B -Vb- —jz — ^ - Qbf np] . (22) 

gsB ^Mb 

Replacing it into Eq. (jH)) and deriving the resulting expression with respect to r we get 
1 /fi 2x2/3 -1/3 dnp qbp 2 , 2 , ^"^^^ 1 9%p dup 

Once again, the extremum in the fermionic density is found by setting dnp/dr = with 
Up ^ 0, which means that the following equation has to be fulfilled: 

-^M^c.> + M^c.,V + ^4 = ' (24) 
gBB Mb 

whose solution again coincides with r = r' as found in the previous Sec. [Eq. (fTHj) ]. 

Due to the attractive boson-fermion interaction, stable trapped ^"K-^'^Rb mixtures may 

only have a limited number of fermions and bosons. If the atom numbers increase above 

n 

some critical values N"^ and Np, an instability occurs|3^]- It has been shown that the mean- 
field framework is . We have calculated 
the stability diagram of the ^°K-^^Rb mixture by solving the coupled mean-field equations 
(jHl) and Q for different values of A''^ and Np. In our study, the instability signature is the 
failure of the numerical iterative process. In particular, the instability onset appears as an 
indefinite growth of the maximum of the densities, which triggers the collapse. 

We display in Fig. 01 the stability diagram for the ^°K-^''Rb mixture in the Nb — Np 
plane. The dots are the theoretical prediction for [Ng,Np). The lines have been drawn 
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to guide the eye and represent the critical instabihty hnes that determine the boundary 
between the stable (left) and unstable (right) regions in four different ) vortex-free 

configurations in concentric trapping potentials (no displacement of the external fermion 
and boson potentials); b) Bose condensate hosting a vortex line with bosons and fermions 
confined by concentric trapping potentials; c) vortex-free configurations with a lO/im dis- 
placement between the trapping potentials; d) Bose condensate hosting a vortex line, plus 
a lOyum displacement between the trapping potentials. 

From Fig. Elone can conclude that, for the ^°K-*^Rb mixture, the presence of a vortex in 
the condensate, or any other mechanism that increases the distance between the maxima of 
the -still- overlapping densities, as for example a sag displacement between the two clouds, 
allows to have a stable mixture for larger particle numbers. The reason is that these mecha- 
nisms diminish the enhancement of the density of both species in the overlap volume caused 
by their attractive mutual interaction. A sag displacement (c) leads to critical numbers 
higher than the presence of a vortex (b), and when both are simultaneously present (d), 
they yield the larger stability region. 

V. SUMMARY 

We have studied Fermi-Bose mixtures with attractive mutual interaction. We have con- 
sidered two potentially interesting cases. In one of them, bosons are in a vortex state, and 
in the other the minima of the trapping potentials for bosons and fermions are shifted from 
each other. 

The position of the maximum of the fermion density is relevant for the collapse of the 
trapped boson-fermion mixture. Indeed, it is around these 'fixed' points that the densities 
increase when the atom numbers approach their critical values, i.e., these are the points 
where the collapse starts. For this reason, simple analytical formulas have been derived to 
describe the effective potential felt by a small fermion amount in the presence of a large 
boson condensate with positive scattering length. These formulas have been used to seek 
the critical points of the fermion density distribution removing the restriction Np << Nb- 
The positions of the critical points only depend on the value of a dimensionless parameter 
(7j) we have introduced in Sec. III. A, whose definition involves the values of the boson and 
fermion scattering lengths, masses and confining frequencies. 
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The validity of these formulas has been assessed comparing the results obtained using 
them with the results obtained from the numerical solution of the mean-field equations (GP 
equation for bosons coupled to the TFW equation for fermions). We have numerically shown 
that the position of the maximum of the fermionic density is very insensitive to the dilution 
value Np/Nsi and have also checked that the analytical values obtained for a low number 
of fermions remain valid up to the critical atomic numbers. 

Finally, we have shown that the critical atomic numbers that the mixture can sustain 
before it collapses may be increased shifting the -still overlapping- confining potentials, or 
when the condensate hosts a quantized vortex line. In the later case, a larger overlap between 
both species has been found, which may favor sympathetic cooling. 
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2l| to derive an expression that 



In this Appendix we use a scaling transformation [IJ, |20 
generalizes the virial theorem and is useful to test the accuracy of the numerical procedure. 
Performing a scaling transformation of the vector position r Ar, it is easy to check that 
to keep the normatization of the order parameter \I'(r ) and of the boson and fermion atomic 
densities, they have to transform as: 



^(f ) ^ ^A(r ) = A3/2^(Af ) 

n{r) nx{r) = \^n{\r) (25) 

Splitting the energy of the system into a kinetic T, a harmonic oscillator Uh, and an inter- 
action part Ug, the above transformations induce the following ones: 
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U,^Ug, = ^g I drnl = X'U, , (26) 

where the last two expressions hold for fermions and bosons as well. The expression for Uh 
supposes a spherically symmetric harmonic trap, the generatization to deformed harmonic 
traps is obvious. 

It is easy to check that in the case of fermions, the kinetic energy in the TFW approxi- 
mation also scales as T — Ta = A^T, and that in the case of boson-fermion mixtures, the 
interaction term Ug^^ scales as Ug^p — > Ug^p^ = X^Uggp. Hence, the total energy of the 
mixture scales as 

E-.E, = X\Ts + Tp) + ^^{Uh, + Uh,) + X\Ug,, + Ugpp + Ug,p) . (27) 

If the scaling is carried out from the equilibrium configuration corresponding to the value 
A = 1, one has 

^ =2{TB + TF)-2{UH, + UH,) + ^Ug,, + Ugpp + Ug,p) = Q . (28) 

A=l 

One sees that at equilibrium, for non-interacting fermions and bosons one recovers the virial 
theorem, namely Tb = and Tp = Uhp- 

When a vortex is present, its superfluid kinetic energy 

Tv 

scales as 

This is quite a natural result in view of the first Eq. ()2fi|l and, as a consequence, Ty may 
be incorporated in the definition of Tb which then represents the total kinetic energy of the 
bosonic component of the mixture. 

If a gravitational sag is considered that displaces the atomic clouds in the z direction, 
its effect can be taken into account changing the confining potential in that direction into 



2M, 



B 



df 



nsir 



n 



(29) 



2M 



B 



df 



nx{r) 



X'Ty 



(30) 
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V{z) = uj'^{z — ZqY. Apart from a trivial constant term proportional to the number of 
fermions and bosons in the trap that is invariant under the scaling transformation, there 
appear new contributions to the total energy of the kind 

Us = -Mzqu^ J drzn{r) (31) 
which scales as Us Us^ = jUs- Eq. then becomes 



^ =2{TB + Tp)-2{UH,+UH,)-iUs,+Us,) + 3iU,,, + U,,, + Ug,,)=0 .(32) 

The above expression is, in a way, a generalization of the virial theorem. We have used 
it to routinely check the accuracy in the numerical solution of the GP and TFW coupled 
equations. Values ~ 0(10) are found for that expression when Uh is ~ O(IO^). This makes 
us confident on the numerical method we have used to find the mean-field structure of the 
atomic mixture. 
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FIG. 1: Boson and fermion density profiles of ^'^K-^^Rb mixtures as a function of z for a cZ^ = 10 /xm 
displacement. The different lines correspond to the profiles for Nb = 10^ and several Np values: 
Np = (dotted line); Np = 10^ (sohd lines); Np = 2.5 x 10^ (dashed lines). The inset shows a 
magnified view of the fermion density distributions. 
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FIG. 2: Boson and fermion density profiles of ^'^K-^^Rb configurations frosting a vortex line with 
m = 1 without displacement of the trapping potentials. In all cases, Nb = 10^. The different lines 
correspond to Np = (dotted line); Np = 10^ (solid lines); Np = 1.5 x 10^ (dashed lines). The 
inset shows a magnified view of the fermion density distributions. 
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FIG. 3: Stability diagram of the ^'^K-^^Rb mixture as a function of the number of atoms. The 
dots are the prediction for the critical number of bosons and fermions. The lines have been drawn 
to guide the eye, and represent the critical instability lines that determine the boundary between 
the stable (left) and unstable (right) regions in four different cases: a) vortex-free configurations 
without displacement of the trapping potentials; b) bosons hosting a vortex line without displace- 
ment of the trapping potentials; c) vortex-free configurations with a 10/xm displacement between 
the trapping potentials; d) bosons hosting a vortex line, plus a 10/xm displacement between the 
trapping potentials. 
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